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We study a model of strongly correlated spinless fermions on a kagome lattice at 1/3 filling, with interactions 
described by an extended Hubbard Hamiltonian. An effective Hamiltonian in the desired strong correlation 
regime is derived, from which the spectral functions are calculated by means of exact diagonalization techniques. 
We present our numerical results with a view to discussion of possible signatures of confinement/deconfinement 
of fractional charges. 
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I. INTRODUCTION 

Excitations with fractional charges have for a number of 
years been a subject of intense investigation in condensed mat- 
ter physics. This strong interest was initially triggered by the 
work of Su, Schrieffer and Heeger—, who showed that in trans- 
polyacetylene such excitations may exist when the system is 
doped^ The fractional charges depend on doping concentra- 
tion and are especially pronounced for certain tt band fillings. 
Their appearance does not require electron-electron interac- 
tions but lattice degrees of freedom play an important role via 
symmetry breaking. Even before the work of Su, Schrieffer 
and Heeger it was realised by Jackiw and Rebbi 3 , that a sys- 
tem with Dirac particles and a symmetry breaking term gives 
rise to excitations with charge e/2. In the mean time the same 
type of theory was applied to graphenei and to a kagome lat- 
tice with 1/3 (2/3) filling A 

A quite different origin of excitations with fractional 
charges was discovered by Laughlin^ with his explanation of 
the fractional quantum Hall effect. Here electron-electron in- 
teractions are crucial. In fact, it is required that the system 
is in the limit of strong electronic correlation where the ki- 
netic energy of the electrons plays a minor role in comparison 
to the Coulomb repulsion. The elimination of the kinetic en- 
ergy is caused by a strong external magnetic field which forces 
the electrons into the lowest Landau level. In distinction to 
the former case, lattice degrees of freedom are not important 
here. It was also found that fractional charges obey fractional 
statistics^ 

Subsequently it was shown that there exists a further class 
of systems with strong electron correlations which yield 
excitations with fractional charges ±e/2. Prerequisites are 
a geometrically frustrated lattice structure like a pyrochlore, 
checkerboard- or kagome lattice as well as certain lattice 
fillings. The excitations with fractional charges can be either 
confined^ or deconfined^ They are low-energy in nature and 
can be treated in the absence of magnetic excitations (for 



instance through the employment of doped dimer models as 
in this paper 1 '). It is interesting that a case of deconfined 
excitations with fractional charges on a three dimensional 
(3D) pyrochlore lattice can be identified. This calls into 
question relationships between fractional charges and frac- 
tional statistics since in 3D there are only fermions or bosons 
possible but not anyons. In contrast, the 2D case readily 
admits the possibility of fractional charges with anyonic 
statistics. Direct comparison of model systems for both 
bosons and fermions 12 - may shed insight on the nature of the 
statistics of emergent fractional charges. 

Emerging experimental techniques lead us to consider also 
how such problems might be investigated in ultracold atomic 
gases; for example, a recent proposal in the form of tun- 
able optical lattice schemes would be readily applicable to 
our models In the context of synthesised materials, spinel 
compounds contain kagome planes and are therefore natural 
candidates for possible experimental verification (a material 
of particular pertinence here is discussed in Ref.fl4l): there is 
furthermore already experimental evidence that electrons in 
pyrochlore lattices can be strongly correlated42>i& 

In the present communication we wish to elaborate on the 
kagome lattice at 1/3 and 2/3 fillings. Since we are interested 
only in charge degrees of freedom, rather than spin degrees 
of freedom, we consider spinless fermions, or their equiva- 
lent, fully spin polarized electrons. In that case a lattice site 
is empty or singly occupied but never doubly occupied due 
to the Pauli exclusion principle. We assume again that the 
electron correlations are strong. We also show that, upon the 
addition or removal of a particle, vacuum fluctuations occur. 
Each of these particles therefore decays into two 'defects', 
each with a fractional charge of ±|. The theory breaks the 
symmetry between the positively and negatively charged ex- 
citations. The paper focuses on the spectral density and its 
properties. In particular we want to show how the strong cor- 
relations and fractional charges show up in that quantity. The 
hole spectral function has symmetry properties which are un- 
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raveled by gauge transformations. It also has features which 
originate from eigenstates of an effective Hamiltonian. The 
latter is obtained in the strongly correlated limit whereby near- 
est neighbour hopping processes are replaced by ring hop- 
ping which lifts the otherwise macroscopic degeneracy of the 
ground state. 

The paper is organized as follows. In Section II, we in- 
troduce the model Hamiltonian. Both the non-interacting and 
interacting cases are considered. An effective Hamiltonian 
valid in the limit of strong electronic correlations is derived 
and found to map directly to a quantum dimer model with res- 
onating plaquettes. — Symmetries of this effective model are 
explained and discussed. In Section EH, the single-particle 
spectral functions are defined and numerical calculations us- 
ing the exact diagonalization method are presented. 

Finally Section IV provides an overview and conclusion. 
Appendix A provides a proof of a generalized sign rule which 
applies in the limit of strong correlations. 



II. MODEL HAMILTONIAN 

We start from an extended Hubbard model on a kagome 
lattice with nearest-neighbor repulsion V. As we are working 
with spinless fermions, i.e., fully spin polarized electrons, no 
on-site repulsion term is necessary. Using second quantized 
notation, the Hamiltonian is written as 

H=-tJ2 ( c h" + Hx -) + V Uin r (1) 
<*.i> 

Here, the operators c 4 (cf ) annihilate (create) a fermion on site 
i. The number operators are given by rii = cjc 4 . The notation 
refers to pairs of nearest neighbors. 

A. Non-interacting case 

If no interactions are present, i.e., V = in Hamiltonian 
Eqn. (Q]), we can simply diagonalize the Hamiltonian in mo- 
mentum space. The spectrum consists of three bands. 5 One of 
them is a flat band ej, = 2t and the other two are dispersive: 
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FIG. 1: (color online) a) The three bands of the kagome lattice are 
shown in the first Brillouin zone, b) The lattice vectors ai and a.2 
are shown on a section of the kagome lattice; they trace out the unit 
cell, of the lattice, c) Three possible configurations that minimise the 
potential energy at 1/3 filling (these constitute just a small fraction of 
possible ground state configurations for this cluster size). 
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FIG. 2: (color online) a) One possible ring exchange process is 
shown here, where particles on sites /, b and c can hop anticlock- 
wise (clockwise) to sites e (d), d (a) and a (e) respectively. Enu- 
merating row-wise as shown, such a process involves always an odd 
number of anticommutations between particles, b) Particle hopping 
around a hexagon as shown always results in a change of the number 
of particles on the sublattice of two; the local gauge transformation 
derivable from this picture explains the symmetry of H e {{. 
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1 ± y/ 4(cos 2 k\ + COS 2 &2 + cos 2 k 3 ) ~ 3 



(2) 

where k n = k • a„ are given with respect to the lattice vectors 
ai = x, &2 = (x + V3y) / 2, = SL2 — sl\ . The lattice vectors 
are defined as depicted in FIG. [TJi. 

FIG. QJ shows the three bands in the first Brillouin zone. 
The three-fold band structure has a lowest energy level with a 
minimum energy of — At. The band picture demonstrates the 
nature of the ground state at 1/3 filling; only the first band is 
filled, so that the Fermi level sits exactly at the Dirac points. 
This property can lead to interesting physics if the Hamilto- 
nian is weakly perturbed as shown in Ref.|5l 



B. Strong interactions 

In this paper we consider mainly the model at 1/3 filling 
in the limit of strong correlation (t/V <C 1). For strongly 
correlated fermion models on certain frustrated lattices it has 
been shown that at low energies fractionalization of charge 
does occur3 In order to study certain effects of that fraction- 
alization, we derive an effective Hamiltonian which acts on 
the low-energy subspace of the full Hilbert space in the strong 
correlation regime. 

A similar technique has been applied to the extended Hub- 
bard model on the checkerboard lattice . 8 ' 19 - 20 
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We first derive an effective Hamiltonian for the exactly 
1/3 filled case ("undoped case") and then consider the case 
in which one particle is added or removed from the system 
("doped case"). Using the effective model, we are able to per- 
form exact diagonalizations of clusters up to 108 sites. 

Undoped case: In order to derive an effective model for 
the quantum strong correlation limit, we start from the limit 
in which t/V = 0. A value of V > will give rise to a local 
constraint, the triangle rule. Configurations \C) which fulfill 
this constraint have exactly one fermion on each triangle; in 
this case the mutual repulsion energy is zero. FIG. [TJ; shows 
different degenerate ground-state configurations. The number 
of degenerate ground states of the system scales exponentially 
as a function of system size. 

As we turn on a small \ t\ <C V, the macroscopic degeneracy 
is lifted. The lowest-order term in t/V that lifts the degener- 
acy stems from ring-hopping processes around hexagons (in- 
deed, ring exchange processes are found in various frustrated 
spin systems^!). All lower order terms lead solely to a con- 
stant energy shift and do not affect the low-energy physics. 
The low-energy excitations, up to third order in t/V, are de- 
scribed by the following effective Hamiltonian 

#e ff = .9WlO}<Ol +H - C -)- (3) 

v J 

This Hamiltonian acts within the manifold of configurations 
which fulfill the triangle rule and the effective ring-hopping 
amplitude is given by g = 12t 3 /V 2 . The growth with system 
size here is also exponential, albeit much slower than that in 
the case of the extended Hubbard model. An approximate 
formulation of this growth is given in 

Ref.lHas (f)* where 
N is the number of occupied sites (particles). An exact value 
for the bulk entropy per particle in the thermodynamic limit 
given in can be compared with the corresponding finite cluster 
quantities.— 

Even though the original Hamiltonian in Eqn. (Q~|) is 
fermionic, the effective Hamiltonian in Eqn. (f3]) has no sign 
change in the matrix elements. Indeed, for all hexagonal flip- 
ping processes it can be shown that 

(^iO><Oi c M- L w 

This remarkable property belonging to all matrix elements 
of the effective Hamiltonian can be obtained through a row- 
wise enumeration of the kagome lattice sites. With such enu- 
meration, we find that an odd number of anticommutations 
must be performed by the c-operators in (3) when evaluating 
the matrix elements, thus always contributing a negative sign. 
This can be seen by considering a ring exchange process on a 
hexagon as labeled in FIG. |2^. In this case, one particle will 
always move to an adjacent site (thereby undergoing no anti- 
commutations which could introduce a sign). The remaining 
two particles in the process both anticommute with a nonzero 
number of particles. Crucially, one must undergo exactly one 
more (or less) anticommutation then the other. Hence, the 
two particle processes always involves an odd number of anti- 
commutations and thus the sign of any such ring exchange 



process is negative, as in Eqn. ©. The argument holds for 
any hexagon on a row-wise enumerated kagome lattice. From 
this observation it follows that the spinless fermion effec- 
tive model at exactly 1/3 filling is equivalent to a hard-core 
bosonic model. Furthermore, the effective model atl / 3 filling 
is equivalent to that at 2/3 filling. 

We remark that the overall sign in Eqn. ([3) is just a matter of 
a simple local gauge transformation. Specifically, the sign of 
g in Eqn. (0) can be changed by multiplying all configurations 
with a phase 

\C)^i^\C), (5) 

where v(C) is the number of particles on the sublattice shown 
in FIG. HJ5. This fact might appear surprising since the actual 
sign of g can typically be gauged away only for the cases of 
bipartite lattices. The sign of g turns out to be inconsequential 
due to the constrained nature of the ring exchange quantum 
dynamics of Eqn. ® and is a desirable feature. A system so 
described by matrix elements which are all non-positive can 
be studied through Monte Carlo simulation. It is also inter- 
esting for physical reasons as it implies that the system must 
have a symmetric eigenspectrum. 



a) b) 




FIG. 3: (color online) a) The honeycomb lattice (black) comprises of 
the links of the medial lattice belonging to the kagome lattice (light 
grey). Occupied sites of the kagome lattice map to dimers on the 
honeycomb lattice, b) Schematic diagram of the dimer representation 
of the so-called "plaquette phase" in which dimers resonate around 
the bold hexagons. 

Configurations |C) on the kagome lattice fulfilling the tri- 
angle rule correspond to dimer configurations on the honey- 
comb lattice (particles are placed here on links, see FIG. [3^. 
Different dimer coverings are orthogonal because any wave- 
function overlap is neglected. In the dimer representation, the 
effective Hamiltonian reads 

^qdm=. 9 e(iO)(Oi+^ c -)- (6) 

v ' 

The phase diagram of the QDM on a honeycomb lattice is to a 
large extent already understood.-^ The case which we consider 
here, i.e., that of Hqbm, has a ground state which corresponds 
to the three-fold degenerate plaquette phase with broken trans- 
lational invariance.^i21 A schematic of this phase is shown in 
FIG. [3}}. Importantly, the separation of a pair of defects in- 
troduced into this phase tends to destroy the plaquette order. 
This, in turn, reduces the number of ring exchange processes 
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FIG. 4: (color online) a) Charge fractionalization due to adding one 
particle to the 1/3 filled lattice, b) Charge fractionalization due to 
removing one particle from the 1/3 filled lattice, c) Each particle or 
hole excitation creates two defects. Each of these lives on one of the 
two sublattices formed by triangles which constitute the kagome lat- 
tice (that is, on a triangular sublattice, like the example shown), d) 
Here one sees two configurations which are connected by a matrix 
element of the hopping operator. The parity of the particle number 
dc(C) on any connecting path C with alternating empty and occu- 
pied sites between the two defects is always even in one case and 
always odd in the other case (two possible paths are shown in each 
case). 



and hence increases the energy of the system. Hence, g acts 
on pairs of defects as a confinement potential; this intriguing 
feature of the plaquette phase results in a ground-state for the 
systems that is both charge ordered and confined. The dy- 
namical nature of the confinement potential is investigated in 
Section[Hl]of this paper. 

Doped case: After the effective low-energy description 
for the undoped case, we now turn to the case where the sys- 
tem is doped, i.e., a particle is either removed (hole-doping) 
or added (particle doping) to a system at 1 /3 filling. 

Placing one additional particle with charge — e onto an 
empty site of a configuration which fulfills the triangle rule 
leads to a violation on two adjacent triangles, see FIG. 0^. 
The energy is increased by 2 V since the added particle has two 
nearest neighbors (charge gap). There is no way to remove the 
violations of the triangle rule by moving particles. However, 
particles on triangles with two particles can hop to another 
neighboring triangle without creating additional violations of 



the triangle rule, i.e., without increase of the repulsion energy 
(see FIG. Hh). By these hopping processes, two local defects 
(=violations of the triangle rule) can separate and the added 
particle with charge — e breaks into two pieces. They carry a 
fractional charge of — e/2 each. In the quantum mechanical 
case (t 0), the separation leads to a gain in the kinetic en- 
ergy of order t. A similar scenario takes places when remov- 
ing one particle from a configuration which fulfills the triangle 
rule. Removing a particles creates two adjacent defects with 
no particle on a triangle. These defects can separate and each 
of them is carrying a charge of +e/2, see FIGHJ}. 

The dynamics of the particle or hole doped system is de- 
scribed by the effective Hamiltonian Eqn. (01 with an added 
hopping term: 

#do P ed= g E(lO>(Ol + H - c - 
v 

t ^P{c\ Cj + H.c.)P. (7) 

Here, the operator P serves to project out high-energy states, 
i.e., it projects out all configurations with additional violations 
of the triangle rule (as we consider only configurations with 
one pair of defects here). 

For the effective Hamiltonian Eqn. ||3}, we showed above 
that the fermionic sign does not influence the low energy ex- 
citations and that the sign of the effective ring-exchange am- 
plitude g is inconsequential. However, as soon as the system is 
doped a fermionic sign problem arises (as the second term in 
(01 does not necessarily preserve the ring-exchange symmetry 
of the first term). It turns out that the projected hopping term 
of defects in the hole doped system (i.e., Hamiltonian Eqn. (0 
with g = 0) is also bipartite. This remarkable feature of the 
hole-doped hopping counters what we would naively expect. 
The defects hop on two triangular sublattices of the honey- 
comb lattice (see FIG. which are not bipartite. However, 
the bipartite nature of this hopping term is in fact a conse- 
quence of the triangle rule to which it is subject; the nature of 
the restricted hopping is clearly distinct from the free hopping 
of particles on the triangular sublattices. Furthermore, the sign 
of t can be changed by simply multiplying all configurations 
\C) with a phase 



\C) 



-1 



\C), 



(8) 



where dc{C) is the number of particles on an arbitrary path C 
between the two defects (see FIG. 2)1), which has alternating 
occupied and empty sites. Note that for finite systems with 
periodic boundary conditions an extra phase also emerges if 
a particle hops around the torus. A detailed derivation of the 
sign rule is given in AppendixlAl 

The projected hopping term as well as the ring exchange 
term for hole defects possesses such a symmetry. However, 
the combination of the two, Hamiltonian i^doped <GK does not 
generally possess a symmetric eigenspectrum (as the ring ex- 
change term symmetry only exists in the absence of the hop- 
ping term). That is, the symmetry no longer holds for the 
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doped system. The reason is that the dynamics of a doped 
hole involves even and odd permutations of fermions, which 
destroy the parity conservation of the sublattice shown in 
FIG. I2J5. Consequently in the presence of hole hopping pro- 
cesses the ring exchange term of the effective Hamiltonian is 
no longer bipartite. 



III. SPECTRAL FUNCTIONS 

The dynamics of the fractional charges can be studied 
through calculation of the spectral function A(\c,uj) = 
A + (k, uS) + ^4~(k, lu)£, for which the particle and hole con- 
tributions are: 



A+(k,u;) 



lim x 



(9) 



k— 7— p ^(<4)W) 

uj + ir\ + tiQ — H 



lim x 

]Tlm«|(c£)t 

a=l 



1 



u> + irj — Eq + H 



(10) 

CklO 



Here, \4>q) is the ground state of a N particle system with a 
corresponding energy E$. The operator c£ annihilates a par- 
ticle with momentum k in band a of the Hamiltonian in the 
non-interacting case and is obtained from the corresponding 
operators in real-space representation by Fourier transforma- 
tion. The spectral functions are directly measured in ARPES 
experiments and are hence a promising candidate for yielding 
experimental signatures of fractionalization. 

Recent studies showed that the quantity A + (k, 10) may also 
be accessible through photoemission spectroscopy undertaken 
in optical lattice experiments.— 

The resulting integrated spectral density is 



(11) 



For a system with a translationally invariant ground state, the 
quantity D(u) can be conveniently calculated in real-space 
representation from local creation and annihilation operators. 

We first consider the DOS for the extended Hubbard model; 
this full treatment is possible as we use our smallest cluster 
with 27 sites. Tuning through the various regimes of the 
model for this quantity, it is possible to follow the transition 
from non-interacting case to the strongly correlated case. 
We can furthermore demonstrate the near equivalence of the 
effective model in the regime of strong correlations (V = 25t 
tunes the model to well within this regime). In FIG. [5] the 
DOS is thus shown for different values of nearest-neighbor 
repulsion V. In the non-interacting case, V = Ot, the spectral 
functions can be calculated directly from the dispersion 



relation Eqn. (0. The contributions result from the allowed 
k-points in the considered cluster. The bandwidth is 6i, in 
agreement with that seen in FIG. Q] 

Moving away from the non-interacting limit, we find that 
a more complex structure emerges. A gap opens between the 
hole and particle contributions; this is due to the 2V energy 
contribution the addition of a particle adds to the system. The 
gap opening signifies a metal-insulator transition and is esti- 
mated to appear at V « 3t£- The delta peaks give way to 
broadened features, characteristic of the high density of low- 
lying excitations that arise in the presence of repulsive inter- 
actions. This is shown for V = lOt. This feature as well as a 
broadened bandwidth, become more marked with increasing 
repulsion energy strength, until we approach the limiting case 
of strong correlations (V is effectively infinite with respect to 
t within this regime and hence all values of V/t here yield 
the qualitative features expected when t/V — > 0). Here the 
broadened structure in the limit of strong correlations has an 
enhanced bandwidth much larger than Qt, indeed it is «12t 
if one considers just the bandwidth of the two contributions 
separately. This can be interpreted as a consequence of the 
fractionalization of the hole or particle defects that occurs in 
the presence of strong repulsive interactions. Furthermore, 
in this limit, the calculated spectral densities for the effec- 
tive model give a valid description of the extended Hubbard 
model, i.e., below the threshold of t/V < 1/25 the effec- 
tive model closely approximates the extended Hubbard model. 
For V = 25t it is seen that the effective and exact (extended 
Hubbard) models closely resemble each other; each peak in 
the exact model has a corresponding feature in the effective 
one. However, with respect to each other they are shown to be 
shifted in energy. This is primarily due to the self-energy con- 
tribution (consisting of all processes where a fermion hops to 
a neighbouring site, or around a triangle, then back to its orig- 
inal site) which is absent in the effective model calculation. 
This discrepancy between the two spectral functions is of or- 
der ~t 2 /V and disappears entirely as t/V — > 0. 



A. Asymptotically free case (g <C t) 

We now start by considering the effective Hamiltonian of 
Eqn. (01 in the limit of infinitesimally small g -C t (i.e., case 
the limit of t/V — > 0). In the effective model, g mediates the 
plaquette ordering strength and therefore acts as a parameter 
for the strength of the confinement of the fractional charges. 
Thus infininitesimally small g can be thought of as the case 
where the fractional charges are effectively deconfined, i.e., 
the spatial extent of the bound pair is effectively infinite 
compared to the cluster size. FIG. |6] shows the k-integrated 
spectral density in this limit for various cluster sizes. As 
cluster size increases, key features of the DOS like bandwidth 
and spectral distribution are maintained. This is expected, as 
the scale of a bound pair is much larger than that of any of 
our clusters. We turn now to the angle-resolved spectra. In 
the following we present results and an analysis pertaining to 
these calculations, for the largest accessible cluster with 108 
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FIG. 5: (color online) The DOS for several ratios of V and t for the extended Hubbard model is shown above. A small cluster with 27 sites 
allows for this full treatment. For ease of comparison, the effective DOS is also plotted for V = 25, i.e., in the strongly correlated regime 
(shaded in grey, it can be seen to closely resemble the full DOS as expected). 
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FIG. 6: (color online) The hole (left panel) and particle (right panel) 
contribution to the local density of states D(uj) for g = 0.01 for 
various cluster sizes. 



possible within our calculations to observe the symmetry 
of the ring exchange term in the undoped case numerically 
(recall from Eqn. (0 that this gauge symmetry arises from the 
bipartiteness of the ring exchange processes). 

The most prominent feature at the T-point in both the hole 
and particle doped spectral functions is a sharp peak (see FIG. 
[7] inset). The peaks originate from the fact shown below that 



IV' 



N-l\ 



=0 |C)and |V iV+1 ) = (4 =0 )W> are ap- 



proximate eigenfunctions of the Hamiltonian Eqn. ([T]i in the 
limit of small |t|/V. The operator c^Lq creates a fermion the 
third band with zero momentum, in the non-interacting case; 
i.e., in the lowest dispersive band, as given by the dispersion 
relations in Eqn. (O. The k-space and real space operators for 
band 3 are related by: 



k=0 




(12) 



sites. The spectrum of the hole spectra is almost perfectly 
symmetric as shown in FIG. |7] This is a numerical manifes- 
tation of the local gauge symmetry, presented in Section ITR 
for the projected hopping of hole defects. As discussed in 
Section HTl i/doped does not generally possess such symmetry, 
but as the spectra in FIG. |7]pertains to the limit of very small 
g, the spectral symmetry is nevertheless preserved to very 
good approximation. Although not presented here it is also 



where c^. =0 annihilates a particle at site j. We now demon- 
strate in detail that li/^ -1 ) is an approximate eigenfunction 
in the hole doped case: 



H\ij> 



N-l\ 



ff4=oiv>o) 

[H,cl =0 ] |^o> + 4 =0 #IV>o) 



(13) 
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FIG. 7: (color online) The Brillouin Zone is shown along with 
all distinct k-points for a 108-site cluster. The particle and hole 
contributions are shown in red (light grey) and blue (dark grey) 
to the left and right respectively along three lines of symmetry. 
The k-points in k-space have the coordinates a = (0,0), b — 
(0,0.19tt),c = (0,0,387r),d = (0,0.587r),e = (0.33tt,0),/ = 
(0.17-7T, 0.487r), g — (0.67-7T, 0), in units of 1/a (where lattice con- 
stant a = |ai/2| = |a.2 /2| , see FIG. [TJ Inset: the band-resolved 
spectral functions are shown for the gamma point; the delta peaks in 
both the hole and particle doped cases correspond to the lowest en- 
ergy band (band 3). The lower curves represent the contributions of 
the other two bands 1 and 2. 

where Eq denotes the ground-state energy of the undoped sys- 
tem. For the kinetic energy contribution to the commutator we 
obtain the following expression: 

3 

[^,»,4=o] = EE[«-io] 

a=l k' 

= 4=o4=o = -4* 4 = o: (14) 

where we have used the dispersion relation Eqn. (f2|. The com- 
mutator involving the repulsion energy is calculated using the 
real space representation: 

[7J„ p ,4=o] = ~7=Y1 ( c i n * + c * n ^ ■ (15) 
(ij) 

In the considered limit, each site has four neighboring sites 
which are occupied with probability t/V due to quan- 
tum fluctuations in \4>o) (see FIG. |8}. Thus the sum over all 
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FIG. 8: (color online) a) Particle removed at 1/3 filling; due to 
strong repulsions V, such a particle is surrounded by unoccupied 
sites. However, due to quantum fluctuations these sites are occupied 
with probability t/V. b) Particle added at 1/3 filling; due to strong 
V repulsions any such particle has exactly two nearest neighbours, 
with the probability t/V of having extra nearest neighbours due to 
quantum fluctuations. 

nearest neighbors of an arbitrary site applied to the ground 
state leads to 

[^,4=0] l^o) = V-jL^A- Cilfo) 

i 

= it^"- 1 }. (16) 

Collecting all terms yields 

H\ii N - 1 ) = {-tt + U + E {) )\^ N - 1 )=ES N ~ 1 )- (17) 

Thus the peak in the spectral function A~(k = 0,ui) is lo- 
cated near ui = (note that only corrections of order t/V are 
taken into account). 

An analogous derivation for the particle doped case shows 
that \i/j n+1 ) is an eigenfunction of H with 

H\t N+1 ) = {el + 2V + 2t + E )\j N+1 ) 
= (2V-2t + E )\i> N+1 ) 

and thus the peak in the spectral function A + (k = 0, u>) is 
located near u = 2V — 2t. These peaks are signatures of 
approximate eigenstates of the full Hamiltonian in the limit of 
small t and large V and are therefore expected to arise in spec- 
tral density results for the full Hamiltonian as well as those for 
the effective Hamiltonian. Calculation of this quantity was 
also performed for the extended Hubbard model for a cluster 
with 27 sites.. Calculations for larger clusters are presently 
not possible because the unrestricted Hilbert space becomes 
unmanageable. 

As previously mentioned, a discrepancy between the effec- 
tive and full models, in the position of the hole and particle 
peaks of order t 2 /V, is expected. 



B. Strongly confined case 

Consider now the behavior of the effective model where g is 
significant compared to t. Crucially, the "plaquette ordering" 
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(co-2V)/t (co-2V)/t (co-2V)/t 



FIG. 9: (color online) Here, the spectral function j4 + (k = 0, oj), is shown over a range of t for clusters with 27 and 48 and 75 sites (shown 
from left to right respectively). It is estimated that at a value of g ~ t the bound pairs are critically confined just within the cluster boundary. 
We interpret the formation of two modes at g ~ 2t as a signal of this critical behaviour (see text for details). 



ring exchange term is now tuned to a value such that the frac- 
tional charges are confined at scales comparable to the cluster 
size. 

In the thermodynamic limit, the confined nature of the 
defects would always be observable in the spectral densities. 
However, the energetic confinement of the fractional charges 
on the lattice is masked on a sufficiently small cluster with 
periodic boundary conditions (this is the case of effectively 
deconfined fractional charges as discussed in the previous 
subsection). For sufficiently large g this is no longer the case; 
we observe bound states at the point at which the confinement 
strength brings the two defects just within the boundaries 
of the cluster. We present numerical results and analyses 
precisely outlining this behaviour. 

As two defects separate from each other through hopping 
processes, they generate a "string" connecting them to one 
another. There is an increase of kinetic energy along the 
string locally, as the separating defects reduce the number 
of flippable hexagons through their movement through the 
plaquette-ordered background. The resulting confinement 



strength between the two particles is quantified in terms of 
the string tension, which is the energy gradient obtained 
as the defects separate. The potential for a single defect is 
given by U (r) = sgr, where s is the string tension, g is the 
ring-hopping strength and r is the distance between the two 
defects. 

An estimate for the discrete energy levels for a bound state 
of two defects, is given by the solution of the ID Schrodinger 
Equation with U(r) as the potential energy term. The solu- 
tions, subject to appropriate boundary conditions, are given 
by: 

E N = -Root Airyjv x (sg) 2 / 3 , (18) 

where RootAiry N is the Ath root of the Airy function Ai{r). 
Comparing these energy levels with values of U(r) used in 
our numerical simulations, we estimate that for g w t, con- 
fined defects in the lowest bound state will be just inside the 
boundaries of our finite clusters. The string tension has been 
estimated through a comparison of ground-state energies for 
clusters with static defects separated at various distances. As 
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expected, estimating the critical value of g using ( fl8T i agrees 
in a good approximation with results in FIG.[9j Here at around 
g « 2t for all three cluster sizes, two distinct modes appear 
for the particle contribution to the spectral function at the V- 
point. We interpret these modes as quantized bound states that 
appear as the confinement strength becomes sufficiently large 
to confine the two particles with the cluster boundaries. 
The appearance of bound states in the spectra at a critical 
value of the potential is in principle observable experimen- 
tally. Indeed, analogously, bound states have been recently 
observed in Ising Chains in experiment^ 

Lastly we comment on the bandwidths seen for the DOS. 
From the spectral densities shown for smallest g in FIG. |6]it 
can be seen that in both cases the bandwidths are wl2i. 

This holds in the aforementioned limit of infinitesimal g; 
however, it can be seen that the bandwidths increase with in- 
creasing g, in accordance with the corresponding relaxation 
of the restricted hopping constraints (this can be observed 
in FIG. [9] most prominently for the smallest cluster with 27 
sites). 



IV. SUMMARY AND OUTLOOK 

We have derived an effective model for spinless fermions 
at 1/3 filling in the limit of strong correlations. This has 
been mapped to a quantum dimer model with resonating 
plaquettes; subsequently a term which dopes the system with 
defects has been added. A local gauge symmetry has been 
found for +e hole hopping. The spectral functions for both 
hole and particle doped systems have been calculated and 
this local gauge symmetry has been observed. The primary 
feature of the DOS has been identified as an approximate 
eigenstate of the full Hamiltonian in the limit of strong 
correlations. The confinement of the fractional charges has 
been confirmed qualitatively through our spectral density 
calculations for large values of g in the effective model. In 
future work we intend to extend our calculations of spectral 
densities to include an analogous model system of hard-core 
bosons at 1/3 filling, with a view to understanding the nature 
of the highly non-trivial statistics of the fractional charges. 
As a realization of our model system using ultracold atoms 
has been proposed, a natural extension of the work would 
be to identify spectral features that would be observable in 
experiment. 
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FIG. 10: (color online) Weighted arrows are used to represent a dimer 
configuration on the honeycomb lattice; dimers are represented by 
bold arrows pointing (with weight 2) and empty links by small arrows 
(weight 1). Arrows of weight 2 (1) point towards A (B), one of the 
two bipartite sublattices. A loop with an odd number of alternating 
occupied and empty links is shown, along with labels for right and 
left turning vertices. 



Appendix A: Generalized sign rule 

In this appendix, we prove the generalized sign rule for hole 
doping in the case of an infinite lattice. We first show that 
every loop with alternating occupied and empty sites on the 
kagome lattice has a length of L = 2k sites with k being an 
odd integer if the "triangle rule" is fulfilled. This fact is then 
used to directly derive the sign rule Eqn. ([8j. 

For convenience, we chose to represent the allowed con- 
figurations on the kagome lattice, i.e., those which fulfill the 
"triangle rule" in the dimer representation on a honeycomb 
lattice (see FIGpSJ)). Since the honeycomb lattice is bipartite, 
we can add orientations (arrows) to the dimers. We replace 
each dimer with an arrow of weight 2 pointing from sublat- 
tice A to B and each empty link with an arrow of weight -1 
pointing from B to A (see FIG. [Tot. In this representation, the 
"triangle rule" implies that the sum of all arrows at each ver- 
tex is zero, in other words, the "flux" through each vertex is 
conserved. Every closed loop with alternating occupied and 
empty kagome sites corresponds to a loop of arrows pointing 
in the same direction with alternating weights 2 and -1, re- 
spectively. As illustrated in FIG. \W[ a loop involves ni left 
turns of —60° and n r right turns of 60° which have to add 
up to 360°. Only at "right turns" an arrow is pointing either 
into or out of the surface surrounded by the loop. Thus flux 
conservation requires that n r is even. Since we consider loops 
with alternating weights on the links, the total number of links 
(= n r + in) has to be even. The three conditions for n r and 
ni can be expressed as 

n r — 7i; = 6, n r = 2/i, n r + ni = 2k (Al) 

with h and k being integer. The above equations imply that 
k = 2h — 3 and thus the length of the loop is L = 2k with k 
being odd. This concludes the proof. 

If we now remove a particle from an allowed configuration, 
the two defects can only hop on loops C with alternating occu- 
pied and empty links. Since the loops contain an odd number 
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of dimers in the undoped case (as shown above), the number 
of particles is now even. Thus, the number of particles which 
are in between the two defects is the same in either direction 
and is changed by one after each hopping process. The sign 
of the hopping term is changed by simply multiplying all con- 
figurations | C) with a phase 

\C)^{-l) d ^\C), (A2) 



where dc(C) is the number of particles in between the two 
defects. If a finite lattice with periodic boundary conditions 
is considered, the transformation needs to take into account a 
phase factor when loops wind around the torus. 
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